EPJ manuscript No. 

(will be inserted by the editor) 



Experimental study of the delayed threshold phenomenon in a 
Class-A VECSEL 

CN . A. El Amili 1 a , G. Gredat 1 , M. Alouini 2 , I. Sagnes , and F. Bretenaker 

O; 

O^l I 1 Laboratoire Aime-Cotton, CNRS-Universite Paris Sud 11, 91405 Orsay Cedex, France 
^ 2 Institut de Physique de Rennes, UMR CNRS-Universite de Rennes I 6251, 35042 Rennes Cedex, France 



<N 



CZ3 

o 



> 
00 
00 

m 



3 Laboratoire de Photonique et Nanostructures, CNRS, 91460 Marcoussis, France 
Received: date / Revised version: date 



Abstract. An experimental study of the delayed threshold phenomenon in a Vertical Extended Cavity 

O 

^ ■ Semiconductor Emitting Laser is carried out. Under modulation of the pump power, the laser intensity 

o 

•i-H 



exhibits a hysteresis behavior in the vicinity of the threshold. The temporal width of this hysteresis is 
measured as a function of the modulation frequency, and is proved to follow the predicted scaling law. 
A model based on the rate equations is derived and used to analyze the experimental observations. A 
frequency variation of the laser around the delayed threshold and induced by the phase-amplitude coupling 
is predicted and estimated. 

1 Introduction losses. The laser threshold can then be considered as the 



| bifurcation point between two distinct domains [3 . Below 

. . The laser threshold is a key concept which characterizes threshold, the zero intensity value is the stable solution 

I a given laser and marks the transition from vanishing in- wnile the non zero intensity value is the unstable solution. 

_ pi tensity towards the lasing regime. Defining the threshold Above threshold, the stabilities of these two solutions are 

is a non trivial issue, which sparks off interesting discus- exchanged. From a more experimental point of view, the 

sions [12. Classically, the threshold is defined as the op- threshold can be considered as the region in which sponta- 

erating point for which the unsaturated gain equals the neous and stimulated emissions have comparable weights. 



a Present address: Institut de Physique de Rennes, UMR Furthermore, this bifurcation point is modified when a 
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problem of the behavior of a laser under gain or loss mod- Indeed, the analysis of the dynamical hysteresis, when the 
ulation around threshold has been theoretically treated control of parameter is modulated, is quite difficult and 
by Mandel and Erneux |4|5|6j . When the laser is swept the measurements are inaccurate in semiconductor lasers, 
from below to above threshold, the laser equations pre- This is due to the very short value of the photon life- 
diet a delayed threshold compared with the first threshold time in such lasers (due to the large intra-cavity losses 
obtained by adiabatically sweeping the laser parameters, and the smallness of the cavity) and also to the signifi- 
A linear sweeping across the bifurcation point also leads cant amplified spontaneous emission which may hide such 
to a modification on the threshold profile itself. Indeed, an hysteresis cycle at low modulation frequency. More- 
the delayed threshold creates a sharp transition from zero over, the spontaneous emission noise causes jitter on the 
intensity to laser oscillation. To the best of our knowl- time at which the laser turn-on. However, recently, Ver- 
edge, this behavior was first noticed in the case of He-Ne tical External Cavity Semiconductor Lasers (VECSELs) 
lasers with swept cavity lengths by Mikhnenko [7 and have been shown to exhibit very long cold cavity decay 
Aronowitz [8] . The first experimental investigations have times together with remarkably low noise levels |16|17|18] . 
been carried out by Sharpf et al. |9] with an Ar + laser Moreover, these lasers can easily operate in the class-A 
in which the intra-cavity losses are linearly swept across dynamical regime [19]. Thus, the aim of this paper is to 
the threshold. A similar effect was studied in two differ- report an experimental study of a VECSEL under modu- 
ent laser dynamics such as in the Class-A lasers |9|13|14j lation of the pump laser power. Indeed, such semiconduc- 
and in the Class-B lasers |10|ll|12|15j . When one of the tor lasers are well known for behaving quite similarly to 
control parameters is modulated, the laser intensity ex- the class-A gas lasers cited above. One can thus expect 
hibits a hysteresis defined by two branches corresponding them to exhibit dynamical hysteresis around threshold, 
to the turning on and turning off the laser. In the vicinity However, these lasers are semiconductor lasers in which 
of the threshold, a time lag between the two branches oc- the phase-amplitude coupling, described by Henry's a pa- 
curs. In the works cited above, it has been shown that the rameter [20J, must play a significant role. This coupling 
delay between the turn-on and turn-off branches depends must thus play a role in the frequency of the light field 
on the modulation frequency according an inverse square in the vicinity of the threshold where the laser does not 
root scaling law |5|9j . reach its steady-state regime. To answer these questions, 

we present in the second section the analytical model de- 

In this article, we present an experimental study of the 

scribing the time evolution of the laser intensity close to 

delayed threshold in a Class-A semiconductor laser. This 

threshold, and we propose a definition of the width of the 

work is a complementary investigation of the one reported 

dynamical hysteresis cycle. In the third section, we corn- 
in a reference [15] on a Class-B semiconductor laser diode. 
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pare the delays extracted from this model with our exper- be approximated by: 
imental measurements. Finally, we report an estimation of dl I f I 

the frequency variation of the laser in the vicinity of the 

delayed threshold. This last equation can be integrated to give the following 



T,~7Z (r(0-T-r(,)-l|. <*> 



solution [5]: 

2 Model e ^(^ +ro-1 ) 

1(0) ^ /satr cav ^ W 

We consider here a single-frequency semiconductor laser in , _ /jA . , 

where V (t) is given by: 

which the gain is produced by quantum wells. The gain is 
linearly swept around the laser threshold while the intra- 



V(t)= fe-t^'^- 1 )^'. (5) 



o 



cavity losses are kept constant. The laser is in the good ^ , ,x , v , , , ^ . , , , ,. , 

J ° V (t) can be linked to the Dawson integral, leading to the 

cavity limit, also called class- A regime, meaning that the £ n , , . . , , 

J ° ' ° following asymptotic expansion when £ 



— » oo: 



photon lifetime r cav is much longer than the carriers life- 
time r c . In this case, the time evolution of the carrier den- u ^ — 2p (pt + q) \ L ^2 (pt + g) 2 ^ 4 (pt + g) 4 
sity can be adiabatically eliminated. The laser rate equa- 



pt(pt+2q) 

v(t) 



(6) 

^ -i i . i r , _ where » = x t^— and a = A / 75— ^ — (rn — 1). The expres- 

tions are then reduced to only one equation 01 evolution V ZTcav V ^ Tcav v ' 

n . 1 1 . , ., sion dU) together with equation © has been used to fit the 

01 the laser intensity: wo M_y 

( x experimental recordings of the evolution of the laser inten- 

1 + -y— ) ' sity. This model is also used to define the delay between 

where r (t) is the pumping parameter which is defined as the "adiabatic" threshold, obtained with slow sweeping, 



the ratio between the unsaturated gain and the losses, and and the delayed threshold. In this paper, we define two 
/ sat is the saturation intensity. The pumping parameter critical times as depicted in Figure H Figure H^a) repro- 
is varied linearly from below (r < 1) to above (r > 1) duces the evolution of the laser intensity versus time in 
laser threshold. Since the pump laser power is linearly vicinity of the threshold. The full line is obtained using 
modulated, we express r (t) as: equation (jl]): the laser intensity sharply passes from zero 

value to the non-zero value. Then, once above threshold, 

r (t) =r + r]t . (2) 

the laser intensity grows linearly with time. The dashed 
The sweep rate 77 is simply given by 2 (r max — 7*0) f m , where line shows the "adiabatic" evolution of the laser when the 
?"o> r max, and f m are respectively the initial pump par am- pumping parameter is slowly swept. This solution crosses 
eter, the final pump parameter, and the modulation fre- the time axis at t* which corresponds to the time at which 
quency. Close to threshold, / <C 7 sa t, and equation (pp) can the laser reaches the "adiabatic" threshold. At t* the losses 
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Fig. 1. (a) Full line: evolution of the laser intensity ver- 
sus time, obtained from equation (j3]) with f m = l kHz and 
r] = 146 s _1 . Dashed line: corresponding "adiabatic" laser 
intensity. It crosses the time axis at the date t*. (b) First 
derivative of the intensity versus time. The maximum de- 
fines T. The delay St is defined as T — t*. 
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Fig. 2. Sketch of the experimental set-up. The VECSEL 
cavity is in a planar-concave configuration. The laser is 
optically pumped with a pump laser diode. The current of 
this latter is modulated with a triangular function signal. 
The laser output is detected with a fast photodiode. The 
photodiode voltage is then recorded by an oscilloscope. 



are equal to the gain [r (t*) = 1], leading to: 



t* 



(1-ro) 
rj 



sis cycles, such as those linked to polarization switching 
in VCSELs under current modulation, or thermal effect 
inducing hysteresis in the laser switch-on and switch-off 
( 7 ) points [2TI22] . 



One can define an another critical time T corresponding 

to the time at which the first derivative of the intensity 3 Experiment 

reaches his maximum as shown in Figure [TJb). The dif- 
ference between these two critical times characterizes the 
width of the dynamical hysteresis: 



3.1 Experimental set-up 



St = T-t*. 



(8) 



The laser used in the present experiment is a single- frequency 
VECSEL which operates at ~1 jam (see Fig. [2]). The planar- 
concave cavity is formed by a 1/2-VCSEL which provides 
According to previous works |5|9|10|ll|14j , the delay the gain and acts as high reflectivity planar mirror, and 
St depends on the modulation frequency, and it can be a spherical output mirror (10 cm radius of curvature and 
scaled as St oc 1/vTm- In the experimental part of the 99% reflectivity). The cavity length is L < 10 cm. In the 
present paper, we compare this prediction with the exper- 1/2-VCSEL, the gain is produced by six InGaAs/GaAsP 
imental measurements performed with a VECSEL. This strained quantum wells (QWs) grown on a Bragg mirror 
will also rule out other mechanisms leading to hystere- which is bonded onto a SiC substrate in order to dissipate 
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the heat towards a Peltier cooler. This multilayered stack the upper one decreases smoothly towards zero, 

is covered by an antireflection coating to prevent from any The value of r cav is obtained by fitting the lower branch of 

coupled cavity effect. The gain is broad (~ 6 THz band- the experimental data using equations (jlj) and (j6]). Using 

width), spectrally flat, and has been optimized to reach the data of Figure [3j we obtain r cav « 60 ns. This value 

a low threshold [T8] . In order to force the VECSEL to is consistent with the length and losses of the cavity. The 

oscillate on a single longitudinal mode, we insert an un- full line in Figure [3] is obtained by numerical integration 

coated glass etalon (200 fim thick) inside the cavity. The of Equation (pQ) with this value. This leads to a very good 

laser is optically pumped at 808 nm. The current supply agreement with the experimental data, 
of the pump laser diode is modulated using a triangular The behavior of the laser intensity in Figure [3] can be 

voltage signal provided by a function generator. Starting physically understood in the following way. For the lower 

below the VECSEL threshold r*o < 1, the pumping rate branch (increasing pumping parameter), the response time 

increases linearly up to the maximum pumping parameter of the laser, which diverges close to threshold due to the 

r max , before decreasing linearly back to its initial value, critical slowing down phenomenon, restrains the laser from 

The optical signal is detected with a fast photodiode (rise- responding fast enough to the fact that it is suddenly 

time ~ 13 ns). The photodiode signal is then recorded with above threshold, leading to a delay in its actual threshold, 

an oscilloscope. On the contrary, for the upper branch (decreasing pump- 
ing parameter), the gain suddenly becomes smaller than 

3.2 Dynamical hysteresis effect the losses while the cavity is filled with light. Thus, the 

intracavity intensity decreases exponentially with a decay 

The Figure [3] reproduces the evolution of the laser power 

time governed by the cavity photon lifetime. 

versus pumping rate. The arrows indicate the direction of 
the evolution. The dots represent the experimental data 

3.3 Width of the dynamical hysteresis cycle 

recorded for a modulation frequency equal to 1 kHz. They 

clearly exhibit the dynamical hysteresis evidenced by the In this section, we present our measurements of the evolu- 

existence of two different branches. The lower branch shows tion of the width St of the hysteresis cycle when the pump 

a sharp increase with a threshold delayed from the dashed power modulation parameters are varied. We remind that 

line which represents the "adiabatic" solution. This lat- we define St as in Figure [TJ According to Mandel [5], Arec- 

ter line crosses the horizontal axis at r = 1 (non-delayed chi et al. [10 , and our numerical simulations of equation 

threshold). The upper branch is the evolution of the laser (pQ), we expect St to scale as the inverse square root of f m . 
intensity when the pumping rate decreases down to its ini- To check this behavior, the laser signal evolution is 

tial value ro ~ 0.98. By contrast with the lower branch, recorded for several values of the modulation frequency 



6 
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Fig. 3. Dots: measurements of the laser signal. The curves 
are the result of the averaging of five successive recordings. 
Full line: theoretical evolution of the laser signal, obtained 
from equation (pQ) with f m = 1 kHz, r] = 140 s _1 , and 
fcav = 60 ns. Dashed line: "adiabatic" solution. 

ranging from 1 kHz to 6 kHz by steps of 500 Hz. For each 
modulation frequency, we record a set of five traces of the 
laser signal, which are averaged. We extract St from this 
average. Figure HJ^a) reproduces the evolution of St versus 
f m . The experimental data are represented by the dots 
and show that the delay decreases when the modulation 
frequency increases. This behavior is consistent with the 
results reported previously for gas lasers |5|10|11] . The full 
line in Figure EJa) corresponds to the theoretical values of 
St obtained by simulating equation (pQ) with r cav = 60 ns. 
We obtain a very good agreement with the experimental 
measurements. Figure HJb) shows the same data on a log- 
log plot. A linear fit leads to a negative slope of — 0.50 ± 
0.01. This is in excellent agreement with the expected St ~ 
1/vTm scaling law. To our knowledge, theses results are 
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Fig. 4. (a) Dots: measured width St of the hysteresis cycle 
as a function the modulation frequency f m . Full line: nu- 
merical simulations, (b) Same data on a log-log plot. Full 
line: linear fit leading to a slope equal to —0.50 ± 0.01. 

the first such measurements carried out with a Class-A 
semiconductor laser. 

3.4 Estimation of the laser frequency chirp at 
threshold 

Up to now, we have seen that our class-A semiconductor 
laser behaves exactly like the gas lasers in which dynami- 
cal hysteresis has been observed up to now. However, one 
salient feature of semiconductor lasers is the existence of 
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a non negligible phase-amplitude coupling effect summa- 
rized in Henry's factor a. Since around the delayed thresh- 
old, the population inversion has a non-linear variation as 
a function of time, we can thus expect the laser phase to 
be affected by this effect. The subsequent frequency vari- 
ation can be estimated by rewriting equation (pQ) for the 
slowly varying complex amplitude A of the laser field: 



dA A . r(t) 

— — = (1 — la) j 

dt 2r cav l V 1 ~\~ ~ r r~ 



1 



(9) 



where we have used the fact that our laser is a Class-A 
laser, where we have chosen the e~ ZUJt time dependence 
for the complex field, and where we suppose for simplic- 
ity that I = |*4 1 2 . The laser frequency deviation 8v with 
respect to the cold cavity eigenfrequency is obtained by 
taking the imaginary part of equation (0, leading to: 

a r (t) 



8v (t) 



47TT cav 1 



I ' 



(10) 



r(t) 



When the laser has reached steady state, one has - 7 

-^sat 

1 leading to the following frequency deviation with respect 
to the cold cavity frequency: 



(11) 



47TT cav 

Figure [5] reproduces the evolutions of the intensity and the 
frequency deviation versus time according to equations. 
(PQ) and ([TP]) , with the same parameters as in Figure [TJ We 
have also taken a = 7. The frequency deviation follows the 
time evolution of the population inversion. It reaches its 
maximum value during the hysteresis cycle, i. e., when the 
population inversion is different from its threshold value, 
i. e., when the laser intensity is very different from its adi- 
abatic value. Later, when the laser intensity is back to its 
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Fig. 5. (a) Time evolution of the laser signal at the mod- 
ulation frequency f m = l kHz. (b) Frequency variation as 
function the time calculated with Henry's parameter a = 
7. 

adiabatic value, the frequency variation reaches its steady 
state value Su st « 21.95 MHz. Then, the population inver- 
sion remains equal to the population inversion at thresh- 
old. Notice that the expected peak-to-peak variation of 
the frequency during the delayed threshold phenomenon 
is of the order of 250 kHz, making it quite difficult to ob- 
serve, particularly because it occurs on short time scales 
(about 10 /jls) and for very low output powers. 

4 Conclusion 

In this paper we have reported an experimental observa- 
tion of the delayed threshold phenomenon in a Class-A 
semiconductor laser. We have shown that as long as the 
delayed threshold phenomenon is concerned, our VEC- 
SEL behaves just like CO2 or Ar + laser systems in which 
this effect had already been studied. For example, we have 
shown experimentally and theoretically that the time lag 
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St scales as l/VTm? where f m is the modulation frequency 11. M. Ciofini, R. Meucci, F.T. Arecchi, Phys. Rev. A 42, 482 

of the pump. Our experimental results are in very good (1990) 

agreement with a model based on the rate equations for 12. S. Balestri, M. Ciofini, R. Meucci, F.T. Arecchi, P. Colet, 

a Class-A laser. We have also predicted the existence of a M - San Mi S ue1 ' S - Balle ' ph y s - Rev - A 44 > 5894 ( 1991 ) 

frequency variation which occurs in the vicinity of the de- 13 ' SH ' Chakmakjian, S. Papademetriou, K. Koch, C.R. 

i a ^ u • + u-wu i +. 4. i • tv.' f Stroud Jr., Phys. Rev. A 40, 1858 (1989) 

layed threshold just while the laser starts lasmg. This ire- ' J ' v J 

. , 14. G. Vemuri, Phys. Rev. A 85, 36 (1991) 

quency variation exists only m semiconductor lasers due 

15. J.R. Tredicce, G.L. Lippi, P. Mandel, B. Charasse, A. 

to the fact a / 0. This phenomenon deserves to be ob- 

Chevalier, B. Picqu, Am. J. Phys. 72, 799 (2004) 

served with a dedicated experimental apparatus. Indeed, 

16. G. Baili, M. Alouini, C. Moronvalle, D. Dolfi, F. Brete- 

according to the simple dependence of 8v versus time with 

naker, I. Sagnes, A. Garnache, Opt. Lett. 32, 650 (2007) 

respect to a, such an apparatus would offer a direct and 

17. G. Baili, F. Bretenaker , M. Alouini, L. Morvan, D. Dolfi, 

accurate way to measure the Henry factor in VCSEL ac- 

I. Sagnes, J. Lightwave Technol. 26, 952 (2008) 

tive media. 

18. A. Laurain, M. Myara, G. Beaudoin, I. Sagnes, A. Gar- 
nache, Opt. Express 17, 9503 (2009) 
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